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ABSTRACT 

We use very precise frequencies of low-degree solar-oscillation modes mea- 
sured from 4752 days of data collected by the Birmingham Solar-Oscillations 
Network (BiSON) to derive seismic information on the solar core. We compare 
these observations to results from a large Monte Carlo simulation of standard 
solar models, and use the results to constrain the mean molecular weight of the 
solar core, and the metallicity of the solar convection zone. We find that only 
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a high value of solar metaUicity is consistent with the seismic observations. We 
can determine the mean molecular weight of the solar core to a very high pre- 
cision, and, dependent on the sequence of Monte Carlo models used, find that 
the average mean molecular weight in the inner 20% by radius of the Sun ranges 
from 0.7209 to 0.7231, with uncertainties of less than 0.5% on each value. Our 
lowest seismic estimate of solar metaUicity is Z = 0.0187 and our highest is 
Z — 0.0239, with uncertainties in the range of 12-19%. Our results indicate that 
the discrepancies between solar models constructed with low metaUicity and the 
helioseismic observations extend to the solar core and thus cannot be attributed 
to deficiencies in the modeling of the solar convection zone. 

Subject headings: Sun: helioseismology - Sun: interior - Sun: abundances 

1. Introduction 

In a series of papers, Asplund et al. (2004, 2005b, 2005c) and AUende-Prieto et al. (2001, 
2002) have revised the spectroscopic determinations of the solar photospheric composition. 
In particular, their results have determined carbon, nitrogen and oxygen abundances to 
be lower by about 25% to 35% than previous determinations (Grcvcssc & Sauval 1998; 
hereafter GS98). The revision of the oxygen abundance leads to a comparable change in 
the abundances of neon and argon. Additionally, Asplund (2000) has also determined a 
somewhat lower value (10%) for the photospheric abundance of silicon compared to the 
GS98 value. As a result, all the elements for which abundances are obtained from meteoritic 
measurements have seen their abundances reduced by a similar amount. These measurements 
have been summarized in Asplund et al. (2005a; hereafter AGS05), and the net result is that 
the ratio of the mass fraction of heavy elements to hydrogen in the Sun is Z/X = 0.0165 
(alternatively, Z = 0.0122), about 28% lower than the previous value, Z/X = 0.0229 {Z = 
0.0169) given by GS98. 

The new low-abundance value for the heavy elements, albeit the result of a much more 
sophisticated modeling of the solar atmosphere, has given rise to discrepancies between he- 
lioseismic observations and predictions from solar models constructed with the low value of 
Z/X. Solar models constructed with the GS98 composition have shown a remarkable agree- 
ment with the solar structure, as determined by helioseismology techniques (Ghristensen- 
Dalsgaard et al. 1996; BahcaU et al. 1997; Morel et al. 1999; Basu et al. 2000). However, 
when the AGS05 composition is adopted in the solar models, the predicted surface helium 
abundance is too low and the convective envelope too shallow. Additionally, the model 
sound-speed and density profiles show a degraded agreement with their solar counterparts 
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when compared to predictions from models that use the older GS98 composition (Bahcall & 
Pinsonneault 2004; Basu & Antia 2004; Bahcall et al. 2005a, 2005b, 2005c, 2006; Delahaye 
& Pinsonneault 2006). 

The discrepancy between the low-Z/X models and the helioseismic observations has led 
to attempts to determine the solar metallicity from helioseismic data alone, just as the solar 
helium abundance was determined by helioseismology. Antia & Basu (2006) used helioseismic 
data to estimate a value for Z in the solar convection zone of 0.0172 ± 0.002, i.e., closer to 
the GS98 value and much larger than the AGS05 value. The uncertainty in their results 
arose from uncertainties in the equation of state, and a lack of data on acoustic modes of 
high angular degree (/ > 200). In this paper we look to the solar core, where uncertainties 
in the physics of the equation of state and opacities are much lower, to try to constrain the 
solar metallicity. 

We make use of solar p-mode data derived from observations made by the ground- 
based Birmingham Solar- Oscillations Network (BiSON; Chaphn et al. 1996). The BiSON 
instruments make disc-averaged observations of the Sun in Doppler velocity. BiSON data 
can be used to determine very precise frequencies of low-Z modes (Z < 3) that can be used to 
probe the solar core. The possibility of using these data to shed light on the solar abundance 
problem was explored by Basu et al. (2007). Basu et al. made very specific combinations 
of the low-/ frequencies, the so-called 'small frequency spacings' and 'frequency separation 
ratios', to compare models with the observations. 

The small frequency spacings of the low-Z p modes are given by the combination. 

dll+2{n) = Vn,l - (1) 

where Vn,i is the frequency of a mode of degree / and radial order n. The fine spacings are 
determined predominantly by the sound-speed gradient in the core. Using the asymptotic 
theory of p modes it can be shown that (see e.g., Christensen-Dalsgaard & Berthomieu 1991) 

where R is the solar radius, and A;(n) is the large frequency spacing given by 

A/ (n) = - (3) 

The large frequency separation depends inversely on the sound-travel time between the 
center and the surface of the Sun. The frequencies Vn,! and 1/^1-1,1+2 are very similar and 
hence are affected in a similar way by near-surface effects. By taking this difference in 
frequency a large part of the effects from the near-surface uncertainties cancels out, making 
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the spacings a useful probe of the deep solar interior and core. Some residual effects do 
nevertheless remain. One way of reducing the effects of the near-surface errors is to use the 
frequency separation ratios. The frequency separation ratios (Roxburgh & Vorontsov 2003; 
Oti Floranes, Christensen-Dalsgaard & Thompson 2005; Roxburgh 2005) are formed from 
the small frequency spacings and large frequency spacings of the modes. The separation 
ratios are then constructed according to: 

Ai{n) Ao(n + l) 

Since both the small and large spacings are affected in a similar manner by near-surface 
effects, these ratios are somewhat independent of the structure of the surface. 

Basu et al. (2007) showed that small spacings and separation ratios for models con- 
structed with the old GS98 composition match the observed BiSON spacings and ratios much 
more closely than do the spacings and ratios of models with the lower AGS05 composition. 
In short, models constructed with higher metallicities compare better with the BiSON data 
than do models constructed with lower metallicities, although the level of agreement dete- 
riorates when the metallicity becomes very large. This indicates that we should be able to 
determine solar metallicity using the spacing and ratio data. 

In this paper, we therefore expand on the work of Basu et al. (2007), and use the small 
spacings and separation ratios from BiSON data to determine the metallicity of the Sun. 
We compare the observed spacings and ratios with spacings and ratios of some 12,000 solar 
models. The models, which were made for an extensive Monte Carlo simulation (Bahcall et al. 
2006), account for all the relevant uncertainties entering standard solar model calculations. 
From this comparison, we show that it is possible to place extremely tight constraints on /ic, 
the mean molecular weight averaged over the inner 20 % by radius (i.e., over most of the solar 
core) of the Sun. Since the mean molecular weight in the core is related to the metallicity 
at the surface (i.e., the convection zone), we can also place reasonably precise constraints on 
the heavy element abundance, Z. Both fic and Z are measures of solar metallicity, albeit for 
different regions of the Sun. 

The rest of the paper is organized as follows. The observed data and the models used 
are described in § [2l In § [3], we compare the BiSON separation ratios with the separation 
ratios of two sequences of models made to test the response of the spacings and ratios to fic 
and Z. In this section we also show how we can obtain a seismic estimate of the solar fi^ 
and Z by comparison of the BiSON and model spacings and ratios. In §|l]we determine the 
systematic errors that arise due to uncertainties in the radiative opacities and the equation 
of state — two quantities are that not amenable to a Monte Carlo type study. In § [5] we 
expand the BiSON-model comparison by analyzing the results for a grand total of 12,000 
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Monte Carlo solar models, computed by Bahcall et al. (2006). This analysis allows us to 
test the impact of changes to several input solar model parameters (in particular changes to 
the mixture of the heavy elements) on the results. Finally, we summarize our results in § Ei 

2. Data and models used 

We have made use of Doppler velocity observations made by the BiSON over the 4752-d 
period beginning 1992 December 31, and ending 2006 January 3. Frequencies were deter- 
mined by fitting resonant peaks in the power spectrum of the complete time series to yield 
estimates of the low-/ frequencies (e.g., see Chaplin et al. 1999). Prior to calculation of the 
small spacings and separation ratios, we removed the solar-cycle shifts from the raw fitted 
low-/ frequencies. Details on the process used to remove the solar-cycle shifts can be found 
in Basu et al. (2007) (see also Chaplin et al. 2005). The BiSON spacings and ratios were 
then constructed from these corrected frequency data, and the uncertainties on individual 
fitted frequencies propagated in the usual manner to give the spacing and ratio uncertain- 
ties. Table [1] contains those solar-cycle-corrected BiSON mode frequencies which were used 
to compute the BiSON frequency spacings and separation ratios analyzed in the paper. 

Basu et al. (2007) showed that the separation ratios depend on the molecular weight, 
but they did not determine the exact dependence of the separation ratios on the average mean 
molecular weight of the core. To do so, we use two very different sets of models. The first set 
of models, which we refer to as the test models, consist of two sequences of ten solar models 
each. One sequence of models was constructed with the relative heavy element abundances 
of GS98, while the second sequence was made with the relative heavy abundances of AGS05. 
To fix the Z/X of a given model in either sequence, the individual relative heavy element 
abundances of GS98 (or AGS05) were multiplied by the same constant factor. This factor 
was then changed from one model to another within the sequence. All models in the two 
sequences were constructed with the same nuclear reaction rates, opacities, equation of state 
and diffusion rates. These models were constructed to test the dependence of the separation 
ratios on the average mean molecular weight, /ic and the total heavy element abundance, 
Z. For reference, we include in this paper tables of low-/ mode frequencies for two standard 
models. Table [2] has frequencies for a model with the exact GS98 abundance; while Table [3] 
has frequencies for a model with the exact AGS05 abundance. 

The second set of models comprised a grand total of 12,000 models created for a Monte 
Carlo study (Bahcall et al. 2006), and we refer to these as the Monte Carlo models. The 
characteristics and methods of computation of the models can be found in Bahcall et al. 
(2006). Here, we summarize the salient points only. 
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For each solar model 19 input parameters were drawn randomly from separate proba- 
bility distributions for every parameter (see Bahcall et al. 2006 for more details). Seven of 
the input parameters were nuclear reaction rates for important low-energy fusion reactions. 
The solar age, luminosity, and diffusion coefficient (rate) were the next three parameters. 
The final nine parameters were the abundances of nine heavy elements: C, N, O, Ne, Mg, Si, 
S, Ar and Fe. For each solar model, radiative opacity tables corresponding to the randomly 
generated composition were computed and used. Radiative opacity tables were generated 
using data and codes provided by the Opacity Project group as presented in Seaton (2005). 
Low-temperature opacities (for temperatures under 10^ K) were from Ferguson et al. (2005). 
Variations of the abundances were dealt with in such a way as to give four sequences of 
Monte Carlo models. Choices had to be made regarding the underlying mixture, and the 
probability distribution for the mixture. Two basic mixtures were used: The GS98 and the 
AGS05 mixtures. Probability distributions were then assigned on the basis of two different 
estimates of the uncertainties in the abundances of the nine individual elements: 'conser- 
vative' (large) uncertainties, based on differences between the abundances of the GS98 and 
AGS05 mixtures; and 'optimistic' (small) uncertainties, based on the uncertainties quoted 
by Asplund et al. (2005a). The content of the four sequences of models may be summarized 
as follows: 

GS-Cons — These 5000 models were made with the 'conservative' (large) abundance 
uncertainties, centered on the GS98 mixture. 

GS-Opt — These 1000 models were made with the 'optimistic' (small) abundance un- 
certainties, centered on the GS98 mixture. 

AGS-Cons — These 1000 models were made with the 'conservative' (large) abundance 
uncertainties, centered on the AGS05 mixture. 

AGS-Opt — These 5000 models were made with the' optimistic' (small) abundance un- 
certainties, centered on the AGS05 mixture. 

We have calculated frequencies of low-degree modes for all the models and constructed 
the small spacings and separation ratios in exactly the same manner as for the observations. 

The model frequencies come from adiabatic calculations. Since non-adiabatic effects arc not 
included, this leads to a well-known mismatch between the absolute values of the adiabatic 
model frequencies and the observed frequencies. However, as noted earlier, differences due 
to these near-surface effects arc reduced significantly by taking frequency differences, and 
using the small spacings and separation ratios. 
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3. Dependence of the separations on metallicity 

In order to parametrize tlie relation between metallicity and the small spacings and 
separation ratios, we begin by comparing the BiSON spacings and ratios with the spacings 
and ratios of the test models. To show how the BiSON-model comparisons were made, 
consider the analysis of the separation ratios. We calculated for each model the differences 
between the observed BiSON ratios, ro2(n) and r 13(11), and the model ratios, 'rQ2(n) and 
r'lsi^), i.e., 

An,,+2(n) = r^+2{n) - r^^^+aW. (5) 

These differences were then averaged over n, for each of the Aro2(?T-) and Ari3(n), to yield 
weighted mean differences, (Aro2) and (Aria): 

(ArM+2) = . (6) 

n 

The formal uncertainties of the BiSON spacings, cr^j ij^jyji)-, were used to weight the averages 
(with the usual uncertainty-squared Gaussian weighting applied). We averaged data over 
the ranges where good determinations of the separation ratios were available, here n = 9 to 
25. 

Fig. [1] shows plots of the weighted mean differences of the separation ratios (in %) versus 
In /ic, the natural logarithm of the average mean molecular weight of the core. The top two 
panels show data for the GS98 set, and the bottom two panels show data for the AGS05 
set. The formal uncertainties on each point, which come from the BiSON data, are not 
plotted. They are 0.046% on each (Aro2); and 0.038% on each (Aria). Fig. [2] shows plots 
of the weighted mean differences versus the natural logarithm of the surface heavy element 
abundances, InZ. Again, the top two panels show data for the GS98 set, and the bottom 
two panels show data for the AGS05 set. 

When a straight line was fitted to the data in each plot (solid lines), the fitting coeffi- 
cients indicated that a linear dependence was a good model for the data. The fit for the Z 
data was in all cases described by: 

(Ar;,;+2) = a/,/+2 + A,i+2 In Z, (7) 

where a;^;+2 and A,/+2 are, respectively, the best-fitting estimates of the intercept and gradient 
of the straight line. A similar straight-line fitting model was used for the /ic data, with In /ic 
used as the independent variable. 
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The quality of the Z fits degraded significantly when Z, rather than In Z, was used as 
the independent variable. Visible departures from a straight line were then observed. The 
fits also deteriorated when the fine spacings, rather than the separation ratios, were used. 
We therefore devote the remainder of the paper to analysis of weighted mean differences 
made from the separation ratios of the BiSON data and the solar models. These weighted 
mean differences have then been used to determine the natural logarithm of Z, as opposed 
to Z itself. We adopted a similar approach to our study of the mean molecular weight in the 
core. Here, we used the In /ic of the models as the independent variable for the plots. Again, 
this was because we found that use of /ic, rather than In/Xc, degraded the quality of the fits 
somewhat, although not as much as in the case of Z. 

With reference to Figures [T] and [H it is not surprising that the two measures of metallicity 
affect the separation ratios in similar ways. The quantities /ic and Z are related, in that 
a higher Z results in a higher /tc. Two models with the same Z, nuclear reaction rates, 
opacities and equation of state can have different values of /ic only if the diffusion rates 
are different in the two models. It is, however, not surprising that the dependence of the 
separation ratios on the two parameters is somewhat different given that /ic also depends 
heavily on the helium abundance in the core. All the models are calibrated to have the 
same radius and luminosity at the solar age, and hence differences in Z generally give rise 
to differences in the core helium abundance. 



3.1. Seismic estimates of solar Z and /ic 

If the observational data are unbiased, the location on the abscissa that marks where 
each best-fitting straight line passes through zero on the ordinate will give us a 'seismic' 
estimate of the average mean molecular weight /ic averaged over the inner 20 % by radius, 
and the surface abundance Z of the Sun. (Note that the Z we refer to is always the present- 
day (solar age) surface Z.) These locations are marked on the various panels of Figures [1] 
and [2] by the intersecting dotted lines. Conclusions drawn from the test models will of course 
neglect any dependence of the differences of the separation ratios on changes to other solar 
model input parameters, including changes to the mixture of the heavy elements. We go 
on to discuss the impact of such changes, and the overall error budget, in later sections. 
Here we show explicitly how a value for, and uncertainty on, Z may be estimated from the 
differences. The same procedures give results for /tc. 

From the best-fitting coefficients, we seek to find InZ where (Ari^/+2) = 0. From 
Equation [71 we therefore have: 

InZ = -a;,/+2M,/+2, 
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So that 

Z = exp[-ai^i+2/Pi,i+2\- (8) 
Application of Equation [H] to the (Aro2) and (Aria) data of each set (GS98 or AGS05) will 
give us four seismic estimates of the solar heavy element abundance. Using GS98 models we 
get Z = 0.01798 and Z = 0.01774 for (Aro2) and (Aris) respectively. The corresponding 
results using the AGS05 models are Z = 0.01617 and Z = 0.01611. A similar analysis 
of fic results in values of 0.7253 and 0.7244 when GS98 models are used, and 0.7260 and 

0. 7255 when AGS05 test models are used. These estimates are listed in the fourth column 
of Tables S] and O along with estimates of the goodness-of-fit and uncertainties. 

We make use of the observed scatter (variance) of the differences (Aro2) and (Aris) 
about their best-fitting straight line to estimate the uncertainty on the seismic estimates 
of solar /ic and Z. In short, we translate the characteristic scatter on the ordinate into an 
implied uncertainty on the abscissa. To determine the uncertainty in Z (the same procedure 
was used for /ic) we first determine the set of residuals about the best-fitting straight line, 

1. e., for each point we compute 

6 {Ar 1^1+2) = (An,i+2) - [a/,/+2 + A,«+2ln2']. 

The variance of these residuals yields an estimate of their la standard deviation, which we 
call cT(Ari ,+2) ■ '^^is characteristic uncertainty on the residuals may be translated into an 
implied uncertainty on InZ via the best-fitting gradient, i.e., 

(T(lnZ) = a(An.,+2>/A,«+2- (9) 

The equivalent la limits on Z are then bounded by exp [InZ + cr(lnZ)] and 
exp [InZ — (T(lnZ')]. Our estimate of the uncertainty is itself uncertain through the uncer- 
tainty on the gradient, A,«+2- Thus we require that N be large enough to ensure that the 
best-fitting gradient (and our 'look-up curve') is well constrained. Here, gradients for the Z 
fits were returned to a fractional precision of better than 1 %. 

The third column of Tables |4] and [5] shows the computed cr^An ,+2) (i^ ^)- '^^^ ^^^^ 
column gives the implied la uncertainties on /i^ and Z, which we call cr(/ic) and a{Z). Here, 
the positive and negative uncertainties were the same, at the level of precision of the data, 
because the fractional uncertainties were so small. The sixth column shows the implied 
precision in the determination of /ic and Z (in %). For each of the GS98 and AGS05 sets, 
we also combined the estimates from (Aro2) and (Aris) - on the assumption the estimates 
are independent - to give the estimates shown in the third and sixth rows of the tables. 

Inspection of the /ic results in Table H] shows that all results (both individual and com- 
bined) are consistent with one another. This is not surprising given the almost direct depen- 
dence of the sound speed, and its derivative, on /ic. The combined estimate obtained with 
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the GS98 models differs from that obtained with the AGS05 models by only la. As we shall 
see in § [5], the uncertainties on the seismic estimates of fic do not increase much when other 
changes to the solar model input parameters are considered. 

Inspection of the Z results in Table[5]shows that the individual, and combined, estimates 
for solar Z are all significantly higher than the 'low' Asplund et al. value of Z ~ 0.0122. 
However, our combined GS estimate (0.01785) and its combined AGS counterpart (0.01611) 
differ from each other by ~ 16cr (combined uncertainty). This difference might at first glance 
be seen as a cause for concern and an indication that systematic errors are much larger than 
the random errors caused by uncertainties in the observed frequencies. However, we go on 
to show in § 5 that other systematic sources of error - arising from the sensitivity of the 
separation ratios to other parameters of the solar models, including the relative mixture of 
the heavy elements - mean that realistic estimates of the uncertainties on Z are actually 
larger in size than the uncertainties given in Table [51 This is in stark contrast to what is 
found for /ic (see previous paragraph). 

Disagreement between the seismic results obtained from the GS98 and AGS98 models 
can to a large extent be understood in terms of the differences in the two mixtures. The 
influence of the relative mixture is clear from the fact that the separation ratios for the 
AGS05 models and the GS98 models are different for the same value of Z; the differences are 
less pronounced for fic- This finding is not difficult to understand. For the calibrated solar 
models used in this work, the dominant contribution to the Z and Hc of each model comes 
from different elements. For Z, the dominant elements, in order of importance, are oxygen, 
carbon, neon and nitrogen. The value of fic is determined by the mass fractions of helium and 
hydrogen in the core. The abundances of hydrogen and helium in the core depend strongly 
on the abundances of heavy metals that contribute to the opacity in the core. These metals, 
again in order of importance, are iron, silicon, sulfur and oxygen. The difference between the 
GS98 and AGS05 mixture lies predominantly in the relative abundances of oxygen, carbon, 
nitrogen and neon, and much less so in the abundances of iron, silicon and sulfur. This 
explains why for the same Z, fic is different for the GS98 and AGS05 models, as can be seen 
from Fig. [3l 

Since the separation ratios depend basically on fi^ and temperature, the location at 
which (Ar;^;_|_2) = occurs at slightly different values for the two sets of models. The slopes 
of the (Ari_i+2)-ln Z curves for the GS98 and AGS05 models (Figures [2]) are also different 
for the same reason. We investigate other sources of systematic errors in § [51 The results 
in § [51 show that realistic estimates of the uncertainties on Z are actually larger in size than 
the uncertainties given in Table [51 



4. Uncertainties due to opacity and equation of state 

The results obtained with the test models have two obvious limitations. First, they do 
not test the impact of changes to the relative mixture of the heavy elements, except to reveal 
that changes in relative abundances matter. And second, the models test the dependence of 
the separation ratios on fi^ and Z for a fixed set of solar model input parameters, with only 
the Z/X varied. The study above does not deal with uncertainties in the seismic estimates 
of solar fic and Z caused by other input parameters, such as nuclear reaction rates, diffusion 
rates, uncertainties in the relative mixtures, etc. We investigate these effects by conducting a 
Monte Carlo study, which is described below in § O This study allowed us to test the impact 
of 19 solar model input parameters on the results. There are, however, two important 
inputs that are not amenable to a Monte Carlo study, and these are radiative opacities and 
the equation of state (EOS). These two quantities cannot be described by a single number 
and hence we are forced to use a different approach to determine the uncertainties in the 
separation ratios caused by uncertainties in opacities and EOS. We use an approach similar 
to that used by Bahcall et al. (2006). 

We determined the uncertainty introduced in (Ar;^/+2) by the opacities as follows. We 
computed a pair of solar models with the same input parameters and EOS, but one was 
made with opacities from the OP project (Badnell et al. 2005) and the other with opacities 
from OPAL (Iglesias & Rogers 1996). For this matched pair of models we get {Ari^i+2)op i 
and {Ari^i^2) OPAL i (here the subscript i denotes the pair of matched models). The unbiased 
estimator for the variance of the difference is 

'/2. (10) 

and we adopt this quantity as the standard deviation squared, af{{/S.ri^i^2) (opacity)). In 
order to obtain a more representative value for cropac((Ar; ^+2)) we averaged this difference 
over a set of = 40 matched pairs of models, where the 19 other input parameters were 
varied for different matched pairs. The final expression we adopt for cropac((A'^«,z+2)) is 



cTopac((An,,+2)) = ^AT-i (opacity)). (11) 

An analogous procedure was used for the EOS, but in this case one model in each pair was 
computed using the 2001 OPAL EOS (Rogers 2001, Rogers & Nayfonov 2002), while the 
other model was computed using the 1996 OPAL EOS (Rogers et al. 1996). 

In addition to the uncertainties on (Aro2) and (Aria) given by Equation [11], we applied 
the same procedure to compute the implied uncertainties for the values of /ic and Z predicted 



i^((Arz,/+2) (opacity)) = (An,/+2)op,i - (ArM+2)oPAL,i 
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by the solar models. Uncertainties in the opacity and EOS will affect the solar model results 
for Z and /ic, i-e. Z and /ic will be somewhat different for the two models in each of the 
matched pairs described above. The results are given in Table O 

The impact on the separation ratios of uncertainty in the radiative opacities is easy 
to understand — changes in opacity cause changes in temperature, which in turn change 
sound speed and its derivative, thereby changing the separation ratios. However, since the 
opacity uncertainties are small in the core, the overall effect is quite modest. The impact of 
uncertainty in the EOS is more important. This might seem surprising, until one realizes 
that the 1996 OPAL EOS did not treat relativistic effects properly at temperatures and 
densities relevant to the solar core. This results in a somewhat deficient core structure (see 
Elliott & Kosovichev, 1998). Since the deficiency is mainly in the core, it will affect the low-/ 
separations used in this work disproportionately and cause larger uncertainties in the seismic 
estimates of solar Z and /ic. The updated 2001 OPAL EOS has the correction put in. Thus 
the uncertainty in the seismic estimates of solar Z and /ic caused by EOS uncertainties may 
be considered to be an upper limit to the EOS effects. 



5. The Monte Carlo study 

In this section we present results, using extended sequences of solar models, which seek 
to address the influence of other input parameters on the seismic solar Z results reported in 
§ 13.11 The bulk of the results come from tests on 12,000 solar models created for a Monte 
Carlo study of the dependence of solar model characteristics on different input parameters 
(see Bahcall et al. 2006). As discussed in §[21 the models have 19 different input parameters 
selected at random from a distribution of the inputs. 

Scatter plots of the weighted mean differences of the separation ratios (in %) versus 
the natural logarithm of /ic are shown in Figures H] and [5] for all four sequences of Monte 
Carlo models. Figures [HI and [7] show the corresponding plots against In Z . Uncertainties 
caused by uncertainties in the radiative opacities and EOS have been included by adding 
to the data random components with normal distributions characterized by the standard 
deviations given in Table [61 The solid lines in each panel are the best-fitting straight lines 
for the data. The dotted lines intersect at the location along each best-fitting line where the 
weighted mean difference is zero. 

We again adopted a linear model for the analysis in the light of: (i) the results on the 
test models (discussed in Section [3] above); and (ii) the observed scatter on the plots, which 
precluded us from imposing a more complicated fitting model. Detailed breakdowns of fitting 
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results for the Monte Carlo sequences are presented in Table [7] (for /ic), and Table [8] (for Z). 
The main results obtained from the Monte Carlo sequences are as follows: 

1. Before turning to discussion of the results on /ic and Z, we consider first the (Aro2) and 
(Aria) of the model sequences. In the case of the AGS-Opt sequence, the distribution 
of (Aro2) is characterized by a mean value of —2.05 % and a standard deviation of 
0.75%. For (Aris) the corresponding values are, respectively, —1.80% and 0.54%. 
These values imply a difference with helioseismology measurements of 2.7a and 3.3a 
for (Aro2) and (Aris) respectively. Note from Table [6]that the EOS has a large impact 
on the total standard deviations of the (Aro2) and (Aris) distributions. As discussed 
previously, we are probably overestimating the EOS uncertainties in the separation 
ratios. Consequently, the 2.7 and 3.3-cr differences should be considered as robust 
upper limits to the real discrepancy. 

On the other hand, for the GS-Opt sequence the mean and standard deviation of 
the (Aro2) distribution are —0.43% and 0.76% respectively, while for (Aria) we get 
—0.32 % and 0.53 % for the mean and the standard deviation. These numbers trans- 
late into differences of only ^ 0.6a with the helioseismology measurements. It is worth 
mentioning that we perform this comparison only for Monte Carlo sequences having 
the optimistic choice of uncertainties, because the aim is to compare helioseismology 
measurements with solar models that adopt compositions (central values and uncer- 
tainties) given by the solar abundance determinations, i.e. GS98 and AGS05. 

2. As in the case of results obtained using the test models, the average mean molecular 
weight in the inner 20 % by radius, fic is determined to much higher precision than Z 
(results in Table [7j). The precision in each of the combined measures of /ic is better 
than 0.5%. Even with the improved precision, the estimates given by analysis of the 
four sequences of models are in excellent agreement with one another and provide 
a very robust determination of Hc- The largest difference between any of the two 
combined measures is significant at only ~ la. The average of our four combined 
seismic estimates is /ic = 0.7226. 

3. From the AGS-Opt sequence of models, i.e. that adopting the AGS05 central values 
and uncertainties for the composition, we find that solar models constructed with the 
AGS05 composition have an average /ic of 0.7088 ±0.0029 (la uncertainty). This value 
differs from our combined seismic estimates of solar /ic, which are given in Table [TJ 
by between 3.1 (for the combined AGS-Cons measure) and 3.7o" (for the combined 
AGS-Opt measure) (where a is determined by adding in quadrature the uncertainties 
from the observations and the uncertainty in the theoretical distribution). On the other 
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hand, in the case of models adopting the GS98 composition, we derive an average value 
of /ic = 0.7203 ± 0.0029 for the optimistic (i.e., lower) uncertainties on the abundance 
inputs. All the combined seismic estimates agree with this value to well within Icr. 

4. All seismic estimates of solar Z are high (see Table [8]). They lie noticeably above the 
low Asplund et al. (2005a) (AGS05) value. This is clear even from visual inspection of 
the plots in Fig. [TJ which shows results for models with heavy-element mixtures based 
on the AGS05 values. The location on the abscissa where the best-fitting lines pass 
through zero lie well to the high-Z side of each cloud of points. None of the AGS-Opt 
models has a mean weighted difference greater than zero. Estimation of Z from these 
data therefore amounts to an extrapolation, rather than the interpolation as is possible 
with the GS-Cons and GS-Opt models. 

From the AGS-Opt sequence of models, which has abundance uncertainties consistent 
with Asplund et al. (2005a), we find that solar models constructed with the AGS05 
composition have an average Z of 0.0125 ± 0.0007 (la uncertainty). This value differs 
from our combined seismic estimates of solar Z, which are given in Table M, by be- 
tween 2.0(7 (for the combined GS-Cons measure) and 4.3cr (for the combined AGS-Opt 
measure) (where a is determined by adding in quadrature the uncertainties from the 
observations and the uncertainty in the theoretical distribution). 

5. All seismic estimates of solar Z are slightly higher than the Grevesse & Sauval (1998) 
(GS98) value, of Z = 0.0169 (see Table [8]). This is clear from visual inspection of 
the plots in Fig. [6l which shows results for models with heavy-element mixtures based 
on the GS values. The best-fitting lines pass through zero in the high-Z parts of the 
distributions of points. Differences between the GS98 value and our combined seismic 
estimates of solar Z lie between 0.5a (for the combined GS-Cons measure) and 2.0a 
(for the combined AGS-Opt measure). 

6. All seismic estimates of solar Z in Table [8] are in good agreement with the seismic 
estimate of Antia & Basu (2006), which was Z = 0.0172 ± 0.0020. The Antia & 
Basu result was quoted with a precision of just over 12%. The precision in our four, 
combined measures of Z (rows 3, 6, 9 and 12 of Table [HD ranges from ~ 12 % to ~ 19 %. 
The largest difference between any of the two combined measures is significant at only 
~ la. 

The observed scatter in Figures S] to [7] deserves discussion, as do the differences between 
the results given by the Monte Carlo sequence of models and the test models from § 121 

The scatter in the computed (Aro2) and (Aris) about the best-fitting straight lines 
- as characterized by the a^^ri 1+2} column 3 of Tables [7] and [8] - far exceeds the sizes 
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of the formal uncertainties on the mean differences found in § |3l Recall that the formal 
uncertainty (a result of uncertainties in the BiSON data) is only 0.046 % on each value of 
(Aro2); and 0.038% on each value of (Aris). The uncertainties on the seismic estimates of 
solar Z and [i^ computed from analysis of the Monte Carlo data are therefore dominated 
by scatter introduced by the input parameter choices for the solar models. As we discuss 
below, the largest contribution to this scatter comes from the relative abundance of heavy 
elements that contribute to Z . It is worth pointing out that while the uncertainties on the 
seismically estimated values of solar Z, obtained with the Monte Carlo sequences, are much 
larger (by over an order of magnitude) than the uncertainties obtained from the test models, 
corresponding differences in the uncertainties on the seismic values of /Xc, i.e., o"(/ic), are just 
a few times larger. 

The Icr standard deviation of the fitting residuals, cr(Ar;;+2>' changes depending on 
whether 'conservative' (large) or 'optimistic' (small) uncertainties are used for the input 
abundances to the solar models. The cr(Ar, ,+2> values for the GS-Cons and AGS-Cons se- 
quences are between ~ 20 and ~ 60 % higher than the corresponding values of the GS-Opt 
and AGS-Opt sequences. This is not surprising — the 'conservative' GS-Cons and AGS- 
Cons sequences cover, by their very nature, a large, 'pessimistic' range of input abundance 
values, thus the scatter due to the uncertainties in the relative abundance is also larger than 
that for the 'optimistic' sets. The uncertainties that arise from using the GS-Cons and the 
AGS-Cons models may therefore be treated as respectable upper bounds for (j{Z) and cr(/ic). 

The gradients and intercepts of the best-fitting lines of the relations between the sep- 
aration ratios and the metallicity of the test models (Figure |2]) are steeper and higher, 
respectively, then those of their Monte-Carlo counterparts (Figures [6] and [7]). The gradients 
are in some cases steeper by as much as a factor of 3. There are also differences between 
one Monte Carlo sequence and another, with the gradients differing by up to ~ 35 %; the 
gradients are constrained typically to much higher precision. However, in spite of these dif- 
ferences, the zero-crossing points, which serve to provide the estimates of solar Z, return 
consistently robust values. 

Again, it is not very difficult to explain why the gradients are higher for the sets of 
test models compared to the sets of Monte Carlo models, and why the gradients differ 
from one Monte Carlo set to another. The answer again lies in the impact of changes 
to the relative mixture of heavy elements. Each set of linear test models has the same 
relative mixture of elements; these models represent the case of perfect correlation between 
all the abundance uncertainties. A given value of Z gives a unique value of /ic, and other 
conditions in the core (e.g., temperature), and thus matches to a unique value of (Aro2) and 
(Aris). The Monte Carlo sequences simulate the opposite situation; here, all the abundance 
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uncertainties are assumed to be independent of one another. The different dependences on 
relative abundances of Z and imply that their variations in the MC sequences become, 
to some degree, uncorrelated. For example, large changes in CNO elements give rise to 
large changes in Z but have a much more modest impact on ^i^ (and consequently on the 
separation ratios). The result is that in the MC sequences the one-to-one relation between 
Z, and other conditions in the core is lost, as is the one-to-one correlation between Z and 
(Aro2) or (Aria). The overall effect is to force the hnear fitting gradients to shallower values 
when a range of mixtures is admitted in the models. There are changes in the details of the 
mixtures between the different Monte Carlo sequences, and hence some (albeit much more 
modest) changes in gradient are seen from one sequence to another. 

It is worth adding that changes to the diffusion rates do not alter significantly the fitting 
gradients. If one takes appropriate subsets of the Monte Carlo sequences, one finds that the 
relation between the rate of diffusion and //c is always the same. And the diffusion rate 
and Z appear to be only very weakly correlated. Changes to the diffusion do not affect 
the relation between the fitting gradients of (Aro2) or (Aris) versus \vlZ or In /ic. The only 
effects seen are changes to the intercept of the fits between the separation ratios and the 
metallicities, which change the seismic estimates of Z, and thus diffusion is a relevant source 
of uncertainty. We find that physical inputs other than the relative mixture, the diffusion 
coefficients, and the EOS have a much smaller effect on the uncertainties of the estimated 
values of solar //c and Z. 



6. Results and Discussion 

We have used the frequencies of low-dcgrcc acoustic oscillations of the Sun, determined 
by the BiSON network over a period of 4752 days, to try and determine the metallicity Z of 
the Sun. We did so by comparing the frequency separarion ratios of a large set of solar models, 
from four different Monte Carlo sequences, with the BiSON observations. Specifically, we 
used a weighted average of the difference of the separation ratios. We find that in addition 
to giving good constraints on the solar metallicity, which by definition is the abundance of 
heavy elements in the solar convection zone, the comparison provides an excellent means to 
determine the mean molecular weight of the solar core. 

The frequency separation ratios of the low-degree acoustic modes are sensitive to the 
conditions in the solar core. By using the weighted mean differences of the separation ratios, 
we have obtained seismic estimates for the mean molecular weight of the solar core (/^c)- 
The seismic estimates are robust and depend only very weakly on the solar models used to 
construct the weighted mean differences. All of the four sequences of Monte Carlo models 
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used in this work allow us to determine /ic to a precision of better than 0.5%, and estimates 
from different sequences are consistent with each other to better than la. This is true even for 
the seismic estimates obtained using solar models with the AGS05 composition. These solar 
models do not agree with hclioseismic results on the solar sound-speed and density profiles, 
surface helium abundance, depth of the convective envelope and frequency separation ratios. 
We have now shown that solar models constructed with the AGS05 composition (central 
values and uncertainties) fail to reproduce the seismically determined average mean molecular 
weight of the solar core by more than 3a. On the other hand, solar models with the older 
GS98 composition have /ic well within la of our seismic estimate of /ic in the Sun. 

We get estimates for the solar metallicity in the range between Z = 0.0187 and Z = 
0.0229 for the four Monte Carlo sequences with uncertainties in the range of 12% to 19% on 
each measurement. All our seismic estimates for metallicity are consistent with each other at 
about the la level, and are higher than the solar metallicity derived by Asplund et al. (2005) 
by between 2.1 and 4.3a. Our estimates arc also consistent with the seismically derived value 
oiZ = 0.0172 ±0.0020 obtained by Antia & Basu (2006). Finally, our estimates are slightly 
higher (by between 0.5 and 2a) than the solar metallicity value recommended by Grevesse 
& Sauval (1998). 

The Antia & Basu (2006) results were obtained by looking at the near-surface ionization 
zones. The signature of interest was the change in the adiabatic index resulting from ion- 
ization of material, which depends on the equation of state and the metallicity. The biggest 
source of systematic error in the results of Antia & Basu was therefore the equation of state. 
Here, we have obtained very similar results for solar Z by looking at a region where very 
different physical inputs matter. The main source of systematic error for this study was 
the relative abundance of the heavy elements. Thus errors in the physical inputs used to 
construct the solar models are not the reason why we obtain consistency high estimates of 
solar Z. As a matter of fact, even when we use \ow-Z models, we still get estimates of solar 
Z that are high. 

When the AGS05 abundances were published, it was soon noticed that solar models 
constructed with those abundances did not agree well with the helioseismically determined 
properties of the Sun. In particular, the most obvious discrepancy was in the position of 
the convection-zone base. It was speculated that this implied that the physical inputs to 
the models, particularly the opacities, were incorrect. The fact that the helium abundance 
in the convection zone was also incorrect lead to the belief that the problems are localized 
to the outer parts of the Sun. Based on the results of this paper, we conclude that this is 
not the case and that the problems extend to the core. It therefore seems unlikely that the 
origin of the discrepancy lies in the simplified modeling of the regions close to the tachocline 
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or in the treatment of convection which is adopted in standard solar model calculations. 
While rotation and its associated mixing are certainly not modeled in the standard solar 
models used in this work, it should be noted that solar models that attempt to account for 
rotation (and some degree of rotationally induced mixing) predict a lower helium abundance 
in the solar core, which translates into a lower average mean molecular weight (Palacios et 
al. 2006). Our results seem to indicate that this would make the disagreement between 
models with the AGS05 composition and hehoseismology even worse. 

This paper utilizes data collected by the Birmingham Solar-Oscillations Network (BI- 
SON), which is funded by the UK Science Technology and Facilities Council (STFC). We 
thank the members of the BiSON team, colleagues at our host institutes, and all others, 
past and present, who have been associated with BiSON. GAV acknowledges the support of 
STFC. SB acknowledges partial support from NSF grant ATM-0348837. AMS is partially 
supported through a John Bahcall Membership and a Ralph E. and Doris M. Hansmann 
Membership at the IAS. 
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Table 1. Solar-cycle corrected BiSON frequencies (in /xHz) used in the paper 
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Table 2. Model frequencies (in /xHz) for exact GS98 abundance 



n 


/ = 





/ = 


1 


/ = 


2 


/ = 


3 


8 






1329, 


.356 


1394, 


.304 


1450. 


.684 


9 


1407, 


.157 


1472, 


.596 


1535, 


.579 


1591, 


.042 


10 


1548, 


.061 


1612, 


.210 


1674, 


.187 


1728, 


.762 


11 


1686, 


.350 


1748, 


.958 


1809, 


.846 


1864, 


.909 


12 


1821, 


.823 


1884, 


.864 


1945, 


.720 


2001, 


.096 


13 


1957, 


.418 


2020, 


.779 


2082, 


.327 


2138, 


.317 


14 


2093, 


.866 


2157, 


.501 


2218, 


.607 


2274, 


.732 


15 


2229, 


.760 


2293, 


.493 


2354, 


.395 


2410, 


.169 


16 


2365, 


.099 


2428, 


.289 


2489, 


.179 


2545, 


.531 


17 


2499, 


.604 


2563, 


.249 


2624, 


.244 


2681, 


.130 


18 


2634, 


.334 


2698, 


.557 


2760, 


.298 


2817, 


.711 


19 


2770, 


.118 


2834, 


.665 


2896, 


.607 


2954, 


.613 


20 


2906, 


.128 


2971, 


.186 


3033 


.263 


3091 


.349 


21 


3042, 


.414 


3107, 


.438 


3169, 


.828 


3228, 


.318 


22 


3178, 


.697 


3243, 


.930 


3306, 


.347 


3365 


.269 


23 


3314, 


.908 


3380, 


.526 


3443, 


.372 


3502, 


.561 


24 


3451, 


.650 


3517, 


.369 


3580, 


.551 


3640, 


.271 


25 


3588, 


.594 


3654, 


.679 










26 


3725, 


.807 















-23- 



Table 3. Model frequencies (in /xHz) for exact AGS05 abundance 
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.618 


2215, 


.457 


2271, 


.133 


15 


2226, 


.755 


2290, 


.126 


2350, 


.930 


2406, 


.628 


16 


2361 


.873 


2425, 


.004 


2485, 


.506 


2541, 


.695 


17 


2496, 


.103 


2559, 


.770 


2620, 


.763 


2677, 


.364 


18 


2631, 


.043 


2695, 


.017 


2756, 


.563 


2813, 


.952 


19 


2766, 


.599 


2831, 


.216 


2892, 


.902 


2950, 


.568 


20 


2902, 


.576 


2967, 


.404 


3029, 


.432 


3087, 


.386 


21 


3038, 


.792 


3103 


.708 


3165, 


.737 


3224, 


.098 


22 


3174, 


.772 


3240, 


.032 


3302, 


.407 


3360 


.995 


23 


3311, 


.132 


3376, 


.466 


3439, 


.208 


3498, 


.371 


24 


3447, 


.692 


3513 


.433 


3576, 


.329 


3635 


.823 


25 


3584, 


.526 


3650, 


.497 










26 


3721, 


.747 
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Table 4. /Xc results for 'linear test' models 



Set 


N 


(%) 






C^(/^c)//^c 
(%) 


Xn-2 


GS-Lin ro2 


10 


0.136 


0.7253 


±0.0012 


0.17 


9.9 


GS-Lin ri3 


10 


0.099 


0.7244 


±0.0010 


0.14 


7.4 


GS-Lin comb. 






0.7248 


±0.0008 


0.11 




AGS-Lin 


10 


0.108 


0.7260 


±0.0009 


0.12 


6.1 


AGS-Lin ri3 


10 


0.089 


0.7255 


±0.0008 


0.11 


6.1 


AGS-Lin comb. 






0.7258 


±0.0006 


0.08 
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Table 5. Z results for 'linear test' models 



Set 


N 


(%) 


Z 


aiZ) 


(%) 


AiV— 2 


GS-Lin ro2 


10 


0.057 


0.01798 


±0.00013 


0.7 


1.7 


pc; T T-,„ 
vjo-J_iin / 13 


1 n 




n ni 77 A 
u.ui / 


_i_n nnnns 
znu.uuuuo 


u.o 


n Q 


GS-Lin comb. 






0.01785 


±0.00007 


0.4 




AGS-Lin ro2 


10 


0.056 


0.01617 


±0.00011 


0.7 


1.7 


AGS-Lin ria 


10 


0.050 


0.01604 


±0.00011 


0.7 


1.9 


AGS-Lin comb. 






0.01611 


±0.00008 


0.5 





Table 6. Effective standard deviations due to uncertainties in the radiative opacities and 

equation of state 



Quantity 


Aro2 


Ari3 


Z 






(%) 


(%) 


(%) 


(%) 


Copac 


0.16 


0.064 


0.083 


0.077 


"^EOS 


0.53 


0.28 


0.047 


0.031 



-26- 



Table 7. //c results for 'Monte Carlo' models 



Set 


N 


^(ArM+2> 
to/ \ 

(%) 








o-(a*c)/a«c 

toy \ 
(%) 


GS-Cons tq2 


5000 


0.82 


0, 


.7234 


±0.0051 


0.7 


GS-Cons ri3 


5000 


0.56 


0, 


.7229 


±0.0033 


0.6 


GS-Cons comb 






0, 


.7231 


±0 0030 


0.4 


GS-ODt Tno 


1000 


0.65 


0, 


.7233 


±0.0049 


0.7 


GS-Opt ri3 


1000 


0.38 


0, 


.7228 


±0.0030 


0.4 


*jk5-wpt comD. 






n 

u, 




ztU.UUzD 




AGS-Cons ro2 


1000 


0.76 


0, 


.7214 


±0.0046 


0.6 


AGS-Cons ri3 


1000 


0.52 


0, 


.7206 


±0.0034 


0.5 


AGS-Cons comb. 






0, 


.7209 


±0.0027 


0.37 


AGS-Opt ro2 


5000 


0.63 


0, 


.7239 


±0.0047 


0.6 


AGS-Opt ri3 


5000 


0.37 


0, 


.7225 


±0.0029 


0.4 


AGS-Opt comb. 






0, 


.7229 


±0.0025 


0.3 
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Table 8. Z results for 'Monte Carlo' models 



Set 






Z 


a{Z) 


a(Z)/Z 


GS-Cons ro2 


5000 


1.03 


0.0189 


+0.0059 

— U.UU40 


28 


GS-Cons ri3 


5000 


0.83 


0.0186 


+0.0052 

— U.UU4tU 


25 


GS-Cons comb. 






0.0187 


+0.0039 
— U.UUoU 


19 


GS-Opt ro2 


1000 


0.72 


0.0197 


+0.0057 


26 


GS-Opt ri3 


1000 


0.52 


0.0192 


+0.0043 

— U.UUoO 


21 


GS-Onf rnmh 






01 Q4 


+0.0035 
-0.0028 


1 6 


AGS-Cons ro2 


1000 


0.90 


0.0233 


+0.0072 
-0.0055 


27 


AGS-Cons ri3 


1000 


0.75 


0.0226 


+0.0062 
-0.0049 


25 


AGS-Cons comb. 






0.0229 


+0.0047 
-0.0037 


18 


AGS-Opt ro2 


5000 


0.72 


0.0220 


+0.0047 
-0.0039 


20 


AGS-Opt ri3 


5000 


0.50 


0.0219 


+0.0038 
-0.0032 


16 


AGS-Opt comb. 






0.0219 


+0.0029 
-0.0025 


12 
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Fig. 1. — Upper two panels: weighted mean differences, (Aro2) (upper left-hand panel) and 
(Aris) (upper right-hand panel), versus In/Xc, plotted for a sequence of 10 solar models. 

All models have the same relative mixture of heavy elements, corresponding the Grevesse 
& Sauval (1998) mixture but differ in the total amount of metals. The solid line in each 
panel is the best-fitting straight line. The dotted lines intersect at the location along each 
best-fitting line where the weighted mean difference is zero. Lower panels: as per upper 
panels, but for a sequence of 10 solar models having a relative mixture of heavy elements 
corresponding to the Asplund et al. (2005) mixture. 
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Fig. 2. — Upper two panels: weighted mean differences, (Aro2) (upper left-hand panel) and 
(Aris) (upper right-hand panel), versus \nZ, plotted for a sequence of 10 solar models. 

All models have the same relative mixture of heavy elements, corresponding the Grevesse 
& Sauval (1998) mixture but differ in the total amount of metals. The solid line in each 
panel is the best-fitting straight line. The dotted lines intersect at the location along each 
best-fitting line where the weighted mean difference is zero. Lower panels: as per upper 
panels, but for a sequence of 10 solar models having a relative mixture of heavy elements 
corresponding to the Asplund et al. (2005) mixture. 
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Fig. 3. — The relation between Z and /ic (the average mean molecular weight in the inner 
20% by radius) for models constucted with the GS98 relative heavy element abundances and 
the AGS05 relative heavy element abundances. 
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Fig. 4. — Upper two panels: weighted mean differences, (Aro2) (upper left-fiand panel) 
and (Aris) (upper right-hand panel), versus In/Xc, plotted for the 5000 Monte Carlo 'GS- 
Cons' solar models. These models were made with the 'conservative' (large) abundance 
uncertainties, centered on the GS98 mixture. The solid line in each panel is the best-fitting 
straight line. The dotted lines intersect at the location along each best-fitting line where the 
weighted mean difference is zero. Lower panels: As per the upper panels, but for the 1000 
Monte Carlo 'GS-Opt' solar models. These models were made with the 'optimistic' (small) 
abundance uncertainties, centered on the GS98 mixture. 
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Fig. 5. — Upper two panels: weighted mean differences, (Aro2) (upper left-liand panel) 
and (Aris) (upper right-hand panel), versus ln//c, plotted for the 1000 Monte Carlo 'AGS- 
Cons' solar models. These models were made with the 'conservative' (large) abundance 
uncertainties, centered on the AGS05 mixture. The solid line in each panel is the best-fitting 
straight line. The dotted lines intersect at the location along each best-fitting line where the 
weighted mean difference is zero. Lower panels: As per the upper panels, but for the 5000 
Monte Carlo 'AGS-Opt' solar models. These models were made with the 'optimistic' (small) 
abundance uncertainties, centered on the AGS05 mixture. Although visually the sohd hues 
on the lower plots do not look hke the best-fitting hues, they are the unbiased hues given by 
a least-squares fit. 
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Fig. 6. — Upper two panels: weighted mean differences, (Aro2) (upper left-fiand panel) 
and (Aris) (upper right-hand panel), versus InZ, plotted for the 5000 Monte Carlo 'GS- 
Cons' solar models. These models were made with the 'conservative' (large) abundance 
uncertainties, centered on the GS98 mixture. The solid line in each panel is the best-fitting 
straight line. The dotted lines intersect at the location along each best-fitting line where the 
weighted mean difference is zero. Lower panels: As per the upper panels, but for the 1000 
Monte Carlo 'GS-Opt' solar models. These models were made with the 'optimistic' (small) 
abundance uncertainties, centered on the GS98 mixture. 
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Fig. 7. — Upper two panels: weighted mean differences, (Aro2) (upper left-liand panel) 
and (Aria) (lower right-hand panel), versus InZ, plotted for the 1000 Monte Carlo 'AGS- 
Cons' solar models. These models were made with the 'conservative' (large) abundance 
uncertainties, centered on the AGS05 mixture. The solid line in each panel is the best-fitting 
straight line. The dotted lines intersect at the location along each best-fitting line where the 
weighted mean difference is zero. Lower panels: As per the upper panels, but for the 5000 
Monte Carlo 'AGS-Opt' solar models. These models were made with the 'optimistic' (small) 
abundance uncertainties, centered on the AGS05 mixture. Although visually the sohd hues 
on the lower plots do not look hke the best-fitting hues, they are the unbiased hues given by 
a least-squares fit. 



